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ABSTRACT 

X 

axial coordinate measured from nozzle throat 

The changes in the nozzle wall surface finish due to 

Y 

distance along heat-flux meter measured from 

roughness may affect the gas to wall heat transfer espe- 
cially if the roughness peaks protrude through the sublayer 


gas-side wall 

of the boundary layer. This paper represents a study on 

y 

distance normal to the wall 

the effects of various degrees of surface roughness on heat 
transfer in two conical converging diverging nozzles with 

y + 

non-dimensional distance normal to the wall 

different convergence angles using heated air. Roughness 
levels up to 325 rms were considered. Roughness causes 

P 

angular position of nozzle instrumentation 

detransition from turbulent flow to take place at a lower 
Reynolds number than that for a smooth wall. In the tur- 

9 

dynamic viscosity 

bulent regime the heat transfer is not noticeably affected 
until the roughness height is greater than a given height 
which seems to correspond to the approximated sublayer 

9 

gas density 

height. 


Subscripts: 

NOMENCLATURE 

ad 

adiabatic wall condition 

A 

nozzle cross-sectional area 

D 

based on diameter 

C f 

local skin friction coefficient 

i 

based on enthalpy 

Cj » Cg, Cg 

constants 

m 

on heat flux meter 

D 

local nozzle diameter 

o 

stagnation condition 

h 

local heat transfer coefficient 

ref 

reference condition 

i 

enthalpy 

s 

static condition 

Pr 

Prandtl number 

t 

throat 

P 

pressure 

w 

wall condition 

q 

heat flux 

oo 

free stream condition 

Re 

Reynolds number 

INTRODUCTION 

rms 

root mean square in micro -inches 


Recent studies have been conducted on the gas -to -wall 



heat transfer of accelerated flows in smooth nozzles with 

T 

total temperature 

different convergence angles. These studies have shown 
that the heat transfer in accelerated flow differs consider- 

t 

static temperature 

ably from that predicted by unaccelerated flow theory and 
experiment. ( 1_4 ) However little consideration was given to 

u 

velocity 

the effects of the wall surface finish on the nozzle heat 
transfer. It is therefore necessary to assess the effect of 

u + 

non-dimensional velocity 

surface roughness on nozzle heat transfer. 
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Acceleration causes a thinning of the boundary layer 
whereas roughness causes an increase in the surface 
height. Therefore the changes in the wall surface finish 
due to roughness may affect the gas -to -wall heat transfer 
especially if the roughness peaks protrude through the sub- 
layer of the boundary layer. Applying this to chemical 
rockets, it is found that after long periods of operation at 
high temperatures the walls of most nozzles become rough- 
er due to erosion and oxidation. These differences in sur- 
face texture introduce questions concerning the applicability 
of the smooth wall heat transfer results to a rocket nozzle. 
Another source of roughness in certain rocket nozzles is a 
ceramic coating which provides insulation between the hot 
propellant-gas and the cooled wall. A knowledge of the 
effects of roughness on the heat transfer is desirable in 
order to determine whether the reduction due to the insula- 
tion is offset by an increase in heat transfer resulting from 
the rougher surface. 

In order to obtain an understanding of the heat transfer 
effects of roughness in accelerated flow one must first con- 
sider independently the behavior of heat transfer in nozzles 
(accelerated flow) with smooth walls. The heat transfer 
studies of Refs. (1) to (4) have shown that the heat transfer 
in accelerated flow differs considerably from that predicted 
by unaccelerated flow theory and experiment. In acceler- 
ated flow there are two distinct regimes of heat transfer 
rate and for both the heat transfer is less than that of con- 
ventional turbulent pipe flow. The high Reynolds number 
regime represents turbulent heat transfer for accelerated 
flows and the low Reynolds number regime suggests a lam- 
inarization phenomenon (sometimes called "detransition”). 
The laminar ization depends on the flow acceleration and is 
further discussed in Ref. (5). 

This paper represents an experimental study on the 
effects of various degrees of surface roughness on heat 
transfer in two converging-diverging nozzles. A 60° half- 
angle of convergence nozzle is compared to a 30° half- 
angle of convergence nozzle for acceleration and roughness 
effects on the gas to wall heat transfer characteristics. 

Both nozzles have a divergence half -angle of 15°. The noz- 
zles were roughened by a sandblasting technique to rough- 
ness levels up to 325 rms. Air at a stagnation temperature 
of 970° R was used as the working fluid. The pressure was 
varied from 30 to 300 psia yielding almost an order of mag- 
nitude range in Reynolds number at each station. The Rey- 
nolds number range at the throat station was 6xl 0 5 to 
5x1 0 6 . These operating conditions made it possible to ob- 
tain heat transfer coefficients and parameters in the tur- 
bulent, transition and laminarization regimes. 


APPARATUS 

The experimental apparatus, shown in Fig. 1, com- 
prised a heat exchanger, diffuser, plenum, pipe inlet, noz- 
zle, and exhaust system. Uncooled- and cooled-wall pipe 
inlets having inside diameters of 6. 5 in. were coupled to 
30° and 60° half-angle of convergence by 15° half-angle of 
divergence water cooled rozzles. These nozzles will 
hereinafter be referred to as simply the 30° and 60° noz- 
zles. 

The adiabatic and cooled inlets have total lengths of 
17.0 and 37.6 in. , respectively. A boundary layer was 
initiated at the leading edge of the inlets by means of a 
bleed flow arrangement. Based on the Reynold s number 


level and surveys in the inlets, the velocity boundary layer 
appeared to become turbulent very close to the leading edge. 
In the cooled inlet, the thermal boundary layer developed 
over a length of 24. 2 in. In tests with the adiabatic inlet, 
the thermal boundary layer started to develop at the nozzle 
entrance. 

The water-cooled 30° and 60° nozzle configurations 
along with tables of the instrumentation sites and coordina- 
tes are presented in Figs. 2 and 3, respectively. Each 
nozzle had a nominal radius of curvature and throat diam- 
eter of 1. 5 in. thus providing a contraction area ratio of 
about 18.8. The expansion ratios for the 30° and 60° noz- 
zles were 25. 3 (Mach 5. 0) and 3. 32 (Mach 2.7), respec- 
tively. 

MEASUREMENTS 

Local heat transfer rates and wall static pressures in 
the nozzles were measured at the stations listed in Figs. 2 
and 3. The heat transfer rates and wall temperatures were 
determined from temperature measurements on Inconel 
plug-type wall heat flux meters which have been described 
in Refs. (2) and (5). Ratios of static -to -total pressure at 
each station have been presented for smooth wall tests in 
Ref. (2). In this study the effect of roughness on the meas- 
ured static pressures was determined to be negligible; 
therefore the pressure ratios presented in Ref. (2) are 
applicable to the present results. 

METHOD OF ROUGHENING 

In this study we are interested in a 1 'natural' * or "uni- 
form” type of roughness where the peaks and valleys are 
connected by a gradual slope rather than a vertical face. 
Nozzle wall roughness values up to 325 rms were obtained 
by a sand or grit blasting technique. The range of rough- 
ness levels was obtained by varying the size and material 
of the abrasive, and the blasting pressure. This method 
left a naturally rough surface with minimal damage to the 
instrumentation. The surface was measured with a rough- 
ness meter giving a root-mean -square value of surface 
height expressed in micro -inches. 

DATA REDUCTION 

The local heat transfer rate q was computed from the 
observed temperature gradient in the heat-flux meters. 

This temperature gradient is described by the Fourier con- 
duction equation which can be integrated over the length of 
the heat meters to give 


-qY = C x 4 + C 2 tm + C 3 (1) 


The constants Ci and C 2 were determined from a 
thermal -conductivity calibration of the heat meter material. 
The values of q and C 3 were obtained from the measured 
temperature on the heat meter, t m , and the corresponding 
thermocouple location Y at two positionson the meter. 

The wall temperature was calculated by setting Y = 0 which 
corresponds to the wall position. 

The errors associated with the heat flux measurements 
for smooth walls have been discussed in Ref. (5). These 
errors were expected to be within ±10 percent of the meas- 


2 



ured heat transfer rates. This error is assumed to apply 
to the present results for rough walls. 

The heat transfer results will be presented in terms of 
a heat transfer coefficient hj and also as a non- 
dimensional grouping St^fPr 0 * 7 . The heat transfer co- 
efficient is given by 


14 =: 9-7— 

J ad “ *w 


( 2 ) 


where the adiabatic enthalpy was calculated from 


of accelerated flow. Both a cooled and an adiabatic inlet 
each having smooth walls, were used in order to determine 
the effects of the thermal history of the fluid on the nozzle 
heat transfer. 

Besides the experimental heat transfer data, Fig. 4 
also shows the standard pipe-flow type of nozzle correla- 
tions. The turbulent correlation 

St ref P r0 - 7 = °.°26 Re'^ (7) 

and the laminar correlation 


i . = i + Pr 1/3 (i - ij 
ad s 'os 


( 3 ) 


The non-dimensional heat transfer grouping St^fPr 0 ■ 7 
is given by 


*ref 


Pr 


.0.7 


-^L_ Pr 0 ' 7 
p ref u » 


( 4 ) 


The Stj-gfPr 0 ' 7 grouping will be presented as a func- 
tion of the Reynolds number given by 


Re D, ref “ 


Pref^D 

Pref 


( 5 ) 


A Prandtl number of 0.71 was assumed in. the above 
equations. The subscript "ref" denotes that properties 
were evaluated at a reference enthalpy condition given by 

i re f = i s + 0 . 5^ - i B ) + 0. 22 Pr 1 ^ . ( 6 ) 


RESULTS AND DISCUSSK)N 

In studying the effects of roughness on heat transfer 
in conical nozzles undergoing different rates of acceleration 
both a 60° and 30° nozzle were used. The results for the 
60° nozzle will be presented in a brief form since they 
were reported fully in Ref. ( 6 ) whereas the results for the 
30° nozzle will be treated more thoroughly. 

Heat Transfer in the 60° Nozzle 


A summary of the heat transfer results for the 60° 
nozzle is found in Fig. 4. Heat transfer in the form of a 
Stanton Prandtl number grouping St re f Pr° • 7 is plotted 
with respect to Reynolds number at the throat station for 
both a cooled and an adiabatic inlet. The advantage of this 
particular non-dimensional presentation is that it shows 
the laminarized and turbulent modes of heat transfer as 
distinct regions which can be identified by a Reynolds num- 
ber in which the reference dimension is the local diameter. 
The diameter is ctosen as a convenient dimension since it 
is approximately proportional to the boundary layer thick- 
ness which, of course, would be a more appropriate di- 
mension if it were known (e.g. , ref. (3)). Therefore it 
must be remembered that two points at the same Reynolds 
number based on diameter, but at different stations, will 
not necessarily have the same heat transfer. The throat 
station is presented because it is the region one is most 
concerned with in a nozzle and the heat transfer is typical 


*ref 


Pr 


0.7 


= 0.29 Re 


-0.5 
D, ref 


( 8 ) 


frame the results very well. Although there now are much 
better methods of predicting nozzle heat transfer these cor- 
relations represent the upper and lower limits of the ex- 
perimental results and will serve as reference levels. How- 
ever, they show no comprehension of acceleration, fluid 
thermal history or surface roughness. 

The smooth nozzle was successively roughened three 
times and the resulting heat transfer parameters using a 
cooled inlet are shown in Fig. 4(a). Two separate regimes 
of heat transfer emerge with a third regime connecting the 
two, as was also observed by the authors of Refs. (3) to (5). 
The mode of heat transfer in the upper regime is turbulent. 
In the lower regime the heat transfer reflects a reduction in 
turbulence commonly called laminarization. Although the 
heat transfer may be close to the values based on the corre- 
lation for a laminar boundary layer there is reason to be- 
lieve the structure of the thermal boundary layer remains 
turbulent for the total plenum pressures used in this re- 
port. The reason being that the temperature profiles 
from Ref. (5) were more turbulent-like than laminar. The 
heat transfer regime connecting the upper and the lower 
corresponds to a transition phenomenon. 

In the turbulence and transition regimes an increase in 
roughness height causes an increase in heat transfer rate. 

At the highest Reynolds number, Rej} t re f» all the experi- 
mental values of heat transfer parameter are in the turbu- 
lent regime. At the first level of roughness the Increase in 
heat transfer is negligible (less than 10 percent), at the 
second level it increases by 17 percent and at the third level 
the increase is 30 percent. As the Reynolds number is de- 
creased turbulent flow undergoes detransition. Roughness 
prolongs the regime of turbulent heat transfer rate and 
causes detransition to take place at a lower Reynolds num- 
ber than it normally would for a smooth wall. At the low 
Reynolds numbers the smooth and the 120 rms wall are in 
the laminarization regime where the effect of roughness on 
heat transfer is negligible. 

In Fig. 4(b) the heat transfer results for the adiabatic 
inlet are presented. They are similar to those of the cooled 
inlet except that for a given Reynolds number the corre- 
sponding value of St^fPr 0 - 7 is higher. The increase in 
St^fPr 0 - 7 is 50 and 25 percent for the turbulent and lam- 
inarized regimes respectively. With the cooled inlet the 
thermal boundary layer begins to grow in the inlet and con- 
tinues growing in the nozzle. With the adiabatic inlet the 
thermal boundary layer does not exist in the inlet and starts 
growing at the entrance of the nozzle. This means that at 
all locations in the rozzle the thermal boundary layer pro- 
duced using the cooled inlet is thicker than that using the 
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adiabatic inlet. (3) The turbulent boundary layer can be 
thought of as a conduction layer which transmits or con- 
ducts heat through it. As a conduction layer, its trans- 
missivity decreases with increasing thickness of the layer. 
Therefore the nozzle heat transfer for tests with the cooled 
inlet is less than the values corresponding to tests with the 
adiabatic inlet. 

In working with pipe flow Nikuradse< 7 ) found that in the 
turbulent regime the effects of roughness become signifi- 
cant when the roughness height is greater than what is often 
referred to as the sublayer height. In Ref. (6), a u + , y+ 
model of the boundary layer profile was used to estimate 
the sublayer hei^it for the 60° nozzle. The sublayer 
height, y, was obtained from the following relation: 

y+ D 

Poo 

y = — — (9) 

f-D 

f x 'ROUGH ^ 

It is recognized that the sublayer height is a quantity 
which is not well defined; however the sublayer is usually 
considered to be the region in which the average viscous 
stress is nearly constant. In Eq. (9) this region was as- 
sumed to be given by y + < 20 . 

In the turbulent regime the heat transfer was not no- 
ticeably affected until the roughness height was in the re- 
gion of or greater than the approximated sublayer height. 
Although good agreement was obtained with the Nikuradse 
hypothesis this conclusion should not be emphasized until 
either actual measurements or very accurate calculations 
are made to find the sublayer height. 

Heat Transfer in the 30° Nozzle 

In Fig. 5, the heat transfer coefficient, hj, is pre- 
sented as a function of axial distance from the nozzle throat 
for a given plenum pressure with the nozzle roughness as 
a parameter. This form of presentation shows the direct 
magnitude change in heat transfer coefficient and its axial 
distribution, which is not clear from the previous figures 
of St re fPr 0,7 vs. Rej) re f. The plenum pressures that 
are used are the same as for the 60° nozzle reported in 
Ref. (6). For the 30° nozzle the range of roughness levels 
was extended to include a value of 55 rms. 

In Fig. 5(a), at a stagnation pressure, p Qf of 300 psia, 
the heat transfer coefficient reaches a maximum in the vi- 
cinity of the throat and then decreases as the flow goes 
supersonic. For a stagnation pressure, p Q , of 75 psia. 

Fig. 5(b) shows the same general axial heat transfer dis- 
tribution except that at the stations just upstream of the 
throat there are large differences in the heat transfer co- 
efficient between the smooth wall and the 180 and 325 rms 
walls. At station 9 the heat transfer of the smooth and 
55 rms walls are in the lower portion of the transition re- 
gime while the heat transfer for the two rough walls are 
turbulent. In Figs. 5(b) and (c) the heat transfer reaches 
a local minimum at station 14 then reverses itself and 
reaches a local maximum at station 15. This same phe- 
nomenon was observed for the 60° nozzle ( 6 ) and also in 
Ref. (4) where the investigator used a smooth nozzle wall 
and varied the stagnation pressure. In all these cases the 


increase in heat transfer occurred where the nozzle exit 
cone meets the throat radius of curvature and this change in 
geometry causes a slight adverse pressure gradient. 

In pipe flow there are cases where the surface may 
seem rough to the senses of sight and feel but behaves as a 
smooth surface. ( 8 ) The rough pipe is then said to be hy- 
draulically smooth. Similarly, we find in Fig. 5 that at 
almost all the nozzle stations the heat transfer coefficient 
for the 55 rms wall is the same as that for a smooth wall. 
Therefore, one can assume that there is a specific rough- 
ness value below which a nozzle can be considered hydrau- 
lically smooth. 

Experimental heat transfer parameters at various sta- 
tions in the 30° and 60° nozzles are presented in Fig. 6. 

The plenum pressure is varied in order to get different 
flow rates which in turn give a range of Reynolds numbers 
at each station. The heat transfer parameter, St^fPr 0 * 7 , 
at each station is presented as a function of Reynolds num- 
ber with the nozzle roughness acting as a parameter. The 
variation of St re fPr 0 *' at stations 9, 11, and 14 for the 
30° nozzle are shown since they represent significant loca- 
tions in the subsonic, sonic and supersonic regimes of 
flow. As expected the heat transfer results fall into re- 
gimes of turbulence, transition and laminarization. The 
heat transfer distributions are very similar to those of the 
60° nozzle of Ref. (6) and will be compared in the following 
section. 


Comparison of the Heat Transfer in the 30° and 60° Nozzles 

In Fig. 6 the heat transfer in the two nozzles are com- 
pared at stations which have nearly the same local nozzle 
diameter. As mentioned in the beginning of this section on 
results and discussion, the local nozzle diameter is chosen 
as the characteristic dimension in the Reynolds number 
rather than the more appropriate dimension of boundary 
layer thickness. Therefore, it must be remembered when 
comparing the results for the 30° and 60° nozzles at a 
given station, that two points at the same Reynolds number 
based on diameter do not in theory have the same Reynolds 
number when boundary layer thickness is used as a char- 
acteristic dimension. However, since the appropriate 
boundary layer thickness Reynolds number is an unknown 
quantity the results are presented in terms of a measurable 
quantity; namely, the Reynolds number based on diameter. 

The increases in the turbulent regime heat transfer 
parameter with respect to the smooth wall for the 60° noz- 
zle are 10, 17, and 30 percent for the respective roughness 
values of 120, 180, and 325 rms. Due to the poor condition 
of some of the inaccessible thermocouples in the 30° noz- 
zle, the increases in heat transfer due to roughness are not 
as consistent as those from the 60° nozzle nor are they con- 
sistent with each other from station to station. For both 
nozzles, roughness causes transition from laminarized 
flow to take place at a lower Reynolds number than that for 
a smooth wall. In the laminarization regime the effect of 
roughness on heat transfer is negligible. Although the 
throat data for the 30° nozzle may not agree with this hy- 
pothesis, the authors feel that this discrepancy is due to an 
instrumentation problem. 

One of the consequences of accelerating flow is the 
phenomenon of laminarization. Therefore by decreasing 
the acceleration, as happens when going from the G0° to 
30° nozzle, the Reynolds number range over which lami- 
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narization takes place should decrease. This is borne out 
in Fig. 6, where the 30° nozzle has fewer points in the 
laminarization region than the 60° nozzle. As a result the 
transition from laminarized flow in the 30° nozzle takes 
place at a lower Reynolds number than in the 60° nozzle. 

For corresponding Reynolds numbers and roughness 
values the heat transfer parameter is greater for the 60° 
nozzle than for the 30°. In the 30° nozzle the axial length 
is greater and consequently the boundary layer is thicker 
at the corresponding diameters. As previously mentioned 
the thermal boundary layer acts as a conduction layer and 
the thicker it is the less heat transfer there is across it. 
Exceptions occur when the 60° nozzle heat transfer is in 
the laminarized regime while corresponding points for the 
30° nozzle are in the transition or turbulent regimes. 

A plot of heat transfer coefficient as a function of local 
nozzle diameter is presented in Fig. 7. Data from both 
nozzles are shown using an adiabatic inlet at a stagnation 
pressure, p Q = 75 psia. For the purpose of clarity, only 
the smooth and 325 rms wall surfaces are shown. The 
heat transfer distribution for the rough surface is typical, 
the 60° nozzle having a higher peak value of h| than the 
30° . Although one would normally expect the same relation- 
ship of heat transfer between the two nozzles for the smooth 
wall as occurred for the rough wall, the results are to the 
contrary. The heat transfer coefficient for the 30° nozzle 
is considerably greater than that for the 60° nozzle at the 
stations in the vicinity of the throat. This is one of those 
exceptional cases that occur when the 60° nozzle heat 
transfer is in the laminarized regime while corresponding 
points for the 30° nozzle are in the transition regime. 

CONCLUSIONS 

An experimental investigation has been performed in 
order to study the effects of roughness on heat transfer in 
conical nozzles. For a roughness range from smooth to 
325 rms the following conclusions can be drawn: 

1. There is a roughness height below which a nozzle 
can be considered to be hydraulically smooth. 

2. Roughness causes detransition from turbulent flow 
to take place at a tower Reynolds number than that for a 
smooth wall. 

3. In the laminarization regime, the effect of rough- 
ness on heat transfer is essentially negligible. 

4. In the turbulent regime the heat transfer does not 
seem to be noticeably affected until the roughness exceeds 
a given height which corresponds to the approximated sub- 
layer height. 
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Figure 1. - Schematic diagram of nozzle heat-transfer facility. 
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STATION 

ANGULAR POSITION, 

P. 

DEG 

AXIAL DISTANCE, 
x, 

IN. 

DIAMETER, 

D, 

IN. 

AREA RATIO, 
A/A t 

PRES- HEAT 
SURE FLUX 
TAP METER 

2 

231 51 

-4. 515 

6.250 

17.548 

3 

129 309 

-3. 515 

5.092 

11.648 

4 

180 

0 

-2.512 

3.934 

6.952 

5 

231 51 

-2. 158 

3.528 

5.591 

6 

283 103 

-1812 

3. 128 

4.395 

7 

334 154 

-1.460 

2. 722 

3.328 

8 

26 206 

-L 110 

2.316 

2.410 

9 

77 257 

-.613 

1.760 

1.392 

10 

129 309 

-.175 

1.510 

1.024 

11 

180 

0 

0 

1.492 

1.000 

12 

231 51 

.130 

1.502 

1.013 

13 

283 103 

.255 

1.540 

1.065 

14 

334 154 

.392 

1.604 

i. 156 

15 

26 206 

.634 

1.732 

1.348 

16 

77 257 

1221 

2.042 

1.873 

17 

180 

0 

2.736 

2.858 

3.669 

18 

1 


5.468 

4.322 

8.391 

19 

1 


8.201 

5.792 

15.070 

20 

♦ 1 


11. 032 

7.320 

24.071 


Figure 2. - Instrumentation for 30° - 15 nozzle. 
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Figure 3. - Instrumentation for 60° - 15° nozzle. 
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TURBULENT 

CORRELATION 


Figure 4. - Experimental heat transfer in the 60° half angle 
of convergence nozzle for various degrees of roughness at 
the throat station. 
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(a) Stagnation pressure, p 0 , 300psia. 
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(b) Stagnation pressure, p 0 , 75 psia. 
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(c) Stagnation pressure, p 0 , 30 psia. 

Figure 5. - Experimental heat transfer distribution in the 30° 
half angle of convergence nozzle for various degrees of rough- 
ness with an adiabatic inlet. 
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(a) Subsonic station. 



(b) Throat station. 



REYNOLDS NUMBER, RE DREF 
(c) Supersonic station. 


Figure 6. - Experimental heat transfer in the 30° 
and 60° half angle of convergence nozzles for 
various degrees of roughness with an adiabatic 
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NOZZLE DIAMETER, D, IN. 


Figure 7. - Effects of roughness on heat transfer coefficient in 
a 30° and 60° half angle of convergence nozzle with an adia- 
batic inlet Stagnation pressure, p 0 , 75 psia. 
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